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We construct a new class of higher dimensional black hole solutions of f(R) theory coupled 
to a nonlinear Maxwell field. In deriving these solutions the traceless property of the energy- 
momentum tensor of the matter filed plays a crucial role. In n-dimensional spacetime the 
energy-momentum tensor of conformally invariant Maxwell field is traceless provided we 
take n = 4p, where p is the power of conformally invariant Maxwell lagrangian. These 
^ . black hole solutions are similar to higher dimensional Reissner-Nordstrom AdS black holes 

but only exist for dimensions which are multiples of four. We calculate the conserved and 
thermodynamic quantities of these black holes and check the validity of the first law of 
black hole thermodynamics by computing a Smarr-type formula for the total mass of the 
solutions. Finally we study the local stability of the solutions and find that there is indeed a 
phase transition for higher dimensional f(R) black holes with conformally invariant Maxwell 
source. 
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I. INTRODUCTION 

A large volume of observational evidences indicate that our universe is now experiencing a 
phase of accelerated expansion. There are two main approaches for explanation of this cosmic 
acceleration. The first approach which try to explain the problem in the framework of general 
relativity, requires the existence of a strange type of energy called "dark energy" whose gravity 
is repulsive and consist an un-clustered component through the universe. The second approach 
is to modify general relativity by adding higher powers of the scalar curvature R, the Riemann 
and Ricci tensors, or their derivatives in the lagrangian formulation. Among the latter attempts 
are Lovelock gravity, braneworld cosmology, scalar-tensor theories like Brans-Dicke one and also 
f(R) theories. A fairly comprehensive review on dark energy models can be seen in f(R) 
theories were proved to be able to mimic the whole cosmological history, from inflation to the 
actual accelerated expansion era 

as 

(see also [4|] for a recent review on f(R) theories). Diverse 
applications of f(R) theories on gravitation and cosmology have been also widely studied 4j|, as well 
as multiple ways to observationally and experimentally distinguish them from general relativity. 

When one considers f(R) theory as a modification of general relativity, it is quite natural to 
ask about black hole existence and its features in this theory. It is expected that some signatures 
of black holes in f(R) theory may differ from the expected physical results in Einstein's gravity. 
Therefore, the investigation on f(R) black holes is of particular interest. Some attempts have 
been done to construct black hole solutions in f(R) theories (see fj-8] and references therein). In 
general, in the presence of a matter field, the field equations of f(R) gravity are complicated and 
it is not easy to find exact analytical solutions. However, if one considers the traceless energy- 
momentum tensor for the matter field, one can extract exact analytical solutions from R + f{R) 
theory coupled to a matter field [9(. Since the energy- momentum tensor of Maxwell and Yang- 
Mills fields are traceless only in four dimensions, therefore black hole solutions from R + f{R) 



theory coupled to the matter field were derived only in four dimensions 



generalized to the four dimensional charged rotating black holes [lfj] , charged rotating black string 



a. 



The studies were also 



111 ] and magnetic string solutions [12j in R + /(i?)-Maxwell theory. However, since the standard 
Maxwell energy-momentum tensor is not traceless in higher dimensions, they failed to derive higher 
dimensional black hole/string solutions from R-\- f{R) gravity coupled to standard Maxwell field. 
A natural question then arises: Is there an extension of Maxwell action in arbitrary dimensions 
that is traceless and hence possesses the conformal invariance? The answer is positive and the 
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conformally invariant Maxwell action was presented as [131 ]. 

S m = - f d n x^d(F^F^) P , (1) 

where p is a positive integer, i.e., p G N. The associated energy- momentum tensor of the above 
conformally invariant Maxwell action is given by 



(2) 



where F = F a pF a/3 is the Maxwell invariant. One can easily check that the above energy momen- 
tum tensor is traceless for n = Ap. The theory of conformally invariant Maxwell field is considerably 
richer than that of the linear standard Maxwell field and in the special case (p = 1) it recovers 
the Maxwell action. It is worthwhile to investigate the effects of exponent p on the behavior of 
the solutions and the laws of black hole mechanics. The motivation is to take advantage of the 
conformal symmetry to construct the analogues of the four dimensional Reissner-Nordstrom (RN) 
solutions, in higher dimensions. Recently, the studies on the black object solutions with a nonlinear 

13-15] and Gauss-Bonnet [l|| gravity have got a lot of attentions. 



Maxwell source in Einstein 



In this work we would like to extend the investigation on the conformally invariant Maxwell 
field to f(R) gravity. We will consider the action (JT]) as the matter source of the field equations 
in R + f(R) theory with constant curvature scalar. Our purpose is to find the analogues of the 
four-dimensional charged black hole solutions of i?+/(i?)-Maxwell theory [9(] in higher dimensional 
spacetime. In contrast to the higher dimensional black holes of Einstein gravity with a conformally 
invariant Maxwell source presented in [lij] which has vanishing scalar curvature R = 0, the space- 
time we construct here in R + f(R) gravity coupled to a nonlinear Maxwell field has a constant 
scalar curvature R = Rq. Our solutions also differ from higher dimensional RN solutions in that 
the electric charge term in the metric coefficient goes as r - (™~ 2 ) while in the standard RN case is 
r — 2(n— 3)_ Also, the electric field in higher dimensions does not depend on n and goes as electric 
field in four dimensions. 

This paper is outlined as follows. In Sec. [IIJ we construct exact spherically symmetric black 
hole solutions of R + f(R) theory coupled to a nonlinear Maxwell field in n = Ap dimensions and 
investigate their properties. In Sec. IIII[ we obtain the conserved and thermodynamic quantities of 
the solutions and verify the validity of the first law of black hole thermodynamics. We also study 
local stability of the solutions in this section. We summarize our results in Sec. IIV1 
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II. FIELD EQUATIONS AND SOLUTIONS 

We consider the action of R + f(R) gravity in n-dimensional spacetime coupled to a conformally 
invariant Maxwell field 

S = [ d n x^~g [R + f(R) - (i^T } , (3) 
JM 

where R is the scalar curvature, f(R) is an arbitrary function of scalar curvature, F^ u = d^A^ — 
d u A^ is the electromagnetic field tensor and is the electromagnetic potential. The equations of 
motion can be obtained by varying action ([3]) with respect to the gravitational field g^ u and the 
gauge field A^, 

iV (i + f'(R)) - \aAR + f(R)) + (a^v 2 - v„v„) f'(R) = r^, (4) 

d^^-gF^F'P' 1 ) = 0, (5) 

where F = F a pF a ^ is the Maxwell invariant and the "prime" denotes differentiation with respect 
to R. In order to obtain the constant curvature black hole solution in f(R) gravity theory coupled 
to a matter field, the trace of stress-energy tensor Tn U should be zero [9(. Hence, two candidates 
for the matter field in four dimensions are Maxwell and Yang-Mills fields. Since the assumption 
of traceless energy- momentum tensor is essential for deriving exact black hole solutions in f(R) 
gravity coupled to the matter field, therefore the solutions exist only for n = 4p dimensions. 
Assuming the constant scalar curvature R = Rq = const., then the trace of Eq. (j4|) yields 

n 

R (1 + f'(Ro)) ~-(Ro + f(Ro)) = 0. (6) 

Solving the above equation for Rq, gives 

_ nf(Ro) _ 2n 
R ° ~ 2f(R ) + 2-n = ^2 Af < °- (7) 

Substituting the above relation into Eq. we obtain the following equation for Ricci tensor 

R^(l + f'(Ro)) - ^Ro (l + f'(Ro)) = T^. (8) 

We are looking for the n-dimensional static spherically symmetric solutions. Motivated by the 
metric of higher dimensional charged black holes in Einstein gravity, we assume the metric has the 
following form 

& 2 = -^(r)^ 2 + ^+r 2 <_ 2 , (9) 
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where df2^_ 2 denotes the metric of an unit (n — 2)-sphere and N(r) is a functions of r which should 
be determined. We are seeking for a purely radial electric solution which means that the only 
non- vanishing component of the Maxwell tensor is F tr . In this case, the Maxwell equations ([5]) can 
be integrated immediately, where, for the spherically symmetric spacetime ([9]), all the components 
of F^ u are zero except F tr : 

F tr = -4=r, (10) 

rp 1p— 1 

where q is an integration constant. Substituting n = 4p in the above relation, the Maxwell field 
becomes 

F tr = 1 (11) 

It is important to note that the electric field in higher dimensions does not depend on n and its 
value coincides with the RN solution in four dimensions. Using metric ([9]) and the Maxwell field 
(jlip , one can show that Eq. ([8]) has a solution of the form 

iV [r) ~ 1 r«~3 + r n-2 X (l + // (jRo )) n (n - 1) r ' [U) 
where m is an integration constant which is related to the mass of the solution. In four dimension 

rj 

(n = 4) the solution recovers the result of [9]. In order to have a real solution we should restrict 
ourself to the dimensions which are multiples of four, i.e., n = 4,8,12,...., which means that p 
should be only positive integer, as we mentioned already. 

Next we study the physical properties of the solutions. To do this, we first look for the curvature 
singularities. A simple calculation shows that the Kretschmann scalar R^ y \ K R^ uXK diverges at 
r = 0, it is finite for r ^ and is proportional to Rq as r — > oo. Therefore, there is a curvature 
singularity located at r = 0. As one can see from Eq. ()12p . the solution is ill-defined for /'(Rq) = 
— 1. Let us consider the cases with f'(Ro) > —1 and f'(Ro) < —1 separately. In the first case where 
f'(Ro) > —1, we have a black hole solution (see Fig. [T]). Indeed, for 1 + f'(Ro) > and Ro < 
black hole can have two inner and outer horizons, an extreme black hole or a naked singularity 
provided the parameters of the solutions are chosen suitably (see Fig. [5]) . In the latter case where 
f'(Ro) < —1, we encounter a cosmological horizon for Rq < 0. Indeed, in this case the signature 
of the spacetime changes and the conserved quantities such as mass become negative, as we will 
see in the next section, thus this is not a physical case and we rule it out from our consideration. 

It is apparent that the spacetime described by solution (fT2j) is asymptotically AdS provided 
we define Rq = —n(n — 1)/^ 2 - However, the solution presented here differ from the standard 



FIG. 1: The function N(r) versus r for m = 2, n = 4, q = 1 and i? = -12. /'(i?o) = -0.5 (bold line) and 
f(Ro) = -2 (dashed line). 
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FIG. 2: The function iV(r) versus r for m = 2, g = 1, n = 4 and i?o 
f'{Ro) = —0.54 (continuous line) and f'(Ro) = —0.8 (dashed line). 



-12. /'(i?o) 



-0.2 (bold line), 



higher dimensional Reissner-Nordstrom AdS (RNAdS) solutions since the electric charge term in 
the metric coefficient goes as r - ( n-2 ) while in the standard RNAdS case is r - 2 ( n - 3 ). i n four 
dimensions, with the following replacement 

■2 



aq 



(H/W) 

R -> 4A 



Q 



(13) 
(14) 



the solution reduces to standard RNAdS black holes for A 



3/1 2 . 



171 ] to calculate the quasilocal mass of 



We use the subtraction method of Brown and York (BY) 
the charged f(R) black hole. Such a procedure causes the resulting physical quantities to depend 
on the choice of reference background. In order to use the BY method one should write the metric 
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in the following form 

dv^ 

ds 2 = -W(r)dt 2 + + r 2 dn 2 n _ 2 . (15) 

Since metric ([9]) has the above form, it is sufficient to choose the background metric to be the 
metric (|15p with 

Wo(r) = V (r) = N (r) = 1 + ^ (16) 

Where we have defined Rq = —n(n — \)jl 2 which show that the solutions are asymptotically AdS as 
we mentioned. It is well-known that the Ricci scalar for AdS spacetime should have this value (see 
e.g. [18]). To compute the conserved mass of the spacetime, we choose a timelike Killing vector 
field £ on the boundary surface B of the spacetime (I15D . Then the quasilocal conserved mass can 
be written as 

M = hj B { ( K <* " Kk ^ - ( K "b - K ° h ab) } ^ ( 17 ) 

where a is the determinant of the metric of the boundary B, K® b is the extrinsic curvature of the 
background metric and n a is the timelike unit normal vector to the boundary B. In the context of 
counterterm method, the limit in which the boundary B becomes infinite (£>oo) is taken, and the 
counterterm prescription ensures that the action and conserved charges are finite. Thus, we obtain 
the mass through the use of the above subtraction method of BY as 

M = (n ~ 8 2 ^- 2 m [1 + f'(Ro)] , (18) 

where Q n -2 is the volume of the unit (n — 2)-sphere. In the limiting case (/'(Rq) = 0), this 
expression for the mass reduces to the mass of the n-dimensional AdS black hole. 

III. THERMODYNAMICS OF CHARGED f(R) BLACK HOLES 

In this section we are going to explore thermodynamics of higher dimensional charged f(R) 
black holes. The Hawking temperature of the black holes can be easily obtained by requiring the 
absence of conical singularity at the horizon in the Euclidean sector of the black hole solutions. 
One obtains the associated temperature with the outer event horizon r = r + as 

1 ( dN(r) \ [1 + f'(R )] (2r 2 (n - 1) + 2l 2 (n - 3)) + (-2g 2 )"/ 4 r 2 -^ 2 

" 4tt V dr J r=r+ nPr[l + f'(R )} ' 1 > 

where we have used equation N(r + ) = for omitting the mass parameter m from temperature 
expression. Next, we calculate the entropy of the black hole. Let us first give a brief discussion 
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regarding the entropy of the black hole in f{R) gravity. To this aim, we follow the arguments 
presented in |19j . If one use the Noether charge method for evaluating the entropy associated with 
black hole solutions in f(R) theory with constant curvature, one finds [fj 



(20) 



where A = A-nr^ is the horizon area. As a result, in f(R) gravity, the entropy does not obey the 
area law and one obtains a modification of the area law. Inspired by the above argument, for the 
n-dimensional charged black hole solutions in R + f(R) gravity, we find the entropy as 
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r n + - 2 £l n -2 



[i + f'(Ro)] ■ 



(21) 



The charge of conformally invariant f(R) black holes can be found by calculating the flux of the 
electric field at infinity, yielding 

n (_2)(™-4)/4g(n-2)/2^ n _ 2 



Q 167^1 + f'(Ro) 

The electric potential <£, measured at infinity with respect to the horizon, is defined by 

where x = @t is the null generator of the horizon. We find 

<*> = + 



(22) 



(23) 



(24) 



Then, we investigate the validity of the first law of thermodynamics for higher dimensional charged 
f(R) black hole. For this purpose, we obtain a Smarr-type formula, namely the mass M as a 
function of extensive quantities 5, and Q. Using the expression for the mass, the entropy and the 
charge given in Eqs. (fl~8j) . (l2~Tj) and ([22]) and the fact that N(r + ) = 0, we find 



M(S,Q) 



n 



45 



-l/(n-2) 



16ttP \\ + f'{Ro) 
' l2+ ( 45 N 2li! - 2 " 



[1 + f'(Ro)] 



45 



1 + f(Ro) 



(25) 



where q is a function of Q according to Eq. (|22|) . One may then regard the parameters 5, and Q as 
a complete set of extensive parameters for the mass M(5, Q) and define the intensive parameters 
conjugate to 5 and Q. These quantities are the temperature and the electric potential 

dS) Q \dQj s 



(26) 
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Numerical calculations show that the intensive quantities calculated by Eq. (|26p coincide with Eqs. 
(I19D and (I24p , Thus, the thermodynamics quantities we obtained in this section satisfy the first 
law of black hole thermodynamics 



In this way we constructed all conserved and thermodynamic quantities of charged f(R) black 
holes and verified the validity of the first law of thermodynamics on the event horizon. 

Finally, we study the local stability of the charged R + f{R) black holes in the presence of 
conformally invariant Maxwell source. In the canonical ensemble, the positivity of the heat capacity 
C Q = T/(d 2 M/dS 2 ) Q and therefore the positivity of (d 2 M/dS 2 ) Q is sufficient to ensure the local 
stability. We have shown the behavior of {d 2 M / 8S 2 )q as a function q and r + for different value of 
n, q and f'(Ro) in figures 3-7 where we have taken Rq = — n(n — l)/l 2 . From figures 3 and 4 we see 
that the system is always thermally stable in four dimensions for different value of the parameters. 
However, in higher dimensions the system has an unstable phase as one can see from Fig. 5. As 
an example, we find that in 8-dimensions black holes are always unstable and there is no phase 
transition (see Fig. 6), while in 12-dimensions the system has a transition from unstable phase to 
stable phase (see Fig. 7). 



In general the field equations of f(R) gravity coupled to a matter field are complicated and it 
is not easy to construct exact analytical solutions. Recently, it was shown [9] that by assuming 



dM = TdS + §dQ. 



(27) 



IV. SUMMARY AND DISCUSSION 
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FIG. 3: The function {d 2 M/dS 2 ) Q versus q for I = 1, n = 4, r+ = 0.8. f(R ) = -0.5 (bold line), f'(R ) = 
(continuous line) and /'(i?o) = 0.5 (dashed line). 
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FIG. 4: The function (d 2 M/dS 2 ) Q versus r+ for I = 1, f(R Q ) = 1 and n = 4. g = 0.4 (bold line), q = 1 
(continuous line) and q = 2 (dashed line). 
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FIG. 5: The function (d 2 M/dS 2 ) Q versus q for Z = 1, r+ = 0.8 and f(R ) = 1. n = 4 (bold line), n = 8 
(continuous line) and n = 12 (dashed line). 
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FIG. 6: The function (d 2 M/dS 2 ) Q versus q for I = 1, n = 8 and r+ = 0.8. .f'(i? ) = -0.5 (bold line), 
/'(-Ro) = (continuous line) and f'(Ro) — 0.5 (dashed line). 
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FIG. 7: The function (d 2 M/dS 2 ) Q versus g for I = 1, n = 12 and r + = 0.8. f'(R ) 
f'{Ro) — (continuous line) and f'(Ro) — 0.5 (dashed line). 
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a traceless energy-momentum tensor for the matter field as well as the constant curvature scalar, 
one can extract some analytical black hole solutions in R + f(R) theory coupled to a matter 
ield. Two examples for the traceless T^ u in four dimensions are Maxwell and Yang-Mills fields 
{J. However, the energy-momentum tensor of Maxwell field is not traceless in higher dimensions. 
Seeking for a traceless energy-momentum tensor in arbitrary dimensions, the authors of [r]] found 
a conformally invariant nonlinear Maxwell action which its energy-momentum tensor is traceless 
in n = 4p dimensions where p is a positive integer. They also studied the black hole solutions in 



Einstein gravity with nonlinear Maxwell field 



1$. 



In this paper we obtained a class of higher dimensional black holes from R + f(R) gravity with 
conformally invariant Maxwell source. The two key assumptions in finding these solutions are: 
(i) the constant scalar curvature R = R® and (ii) the traceless energy momentum tensor. These 
solutions are similar to higher dimensional RNAdS black holes with appropriate replacement of 
the parameters, but only exist in dimensions which are multiples of four. Besides, the solutions 
presented here differ from higher dimensional RNAdS black holes in two features. First, the electric 
charge term in the metric coefficient goes as r~( n-2 ) while in the standard RNAdS case is r ~ 2 ( n ~ 3 ). 
Second, the electric field in higher dimensions does not depend on n and goes as electric field 
in four dimensional RNAdS black holes. Our solutions also differ from the higher dimensional 
black holes of Einstein gravity with a conformally invariant Maxwell source [r3] in that they have 
vanishing scalar curvature R = 0, while the obtained solutions here in R + f(R) gravity coupled 
to a nonlinear Maxwell field have a constant curvature scalar R = Rq. In addition, the conserved 
and thermodynamic quantities computed here depend on function f'(Ro) and differ completely 
from those of Einstein theory in AdS spaces. Clearly the presence of the general function f'(Ro) 
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changes the physical values of conserved and thermodynamic quantities. Furthermore, unlike 
Einstein gravity, for the black hole solutions obtained here in f(R) gravity, the entropy does not 
obey the area law. 

After studying the physical properties of the solutions, we computed the mass, charge, electric 
potential and temperature of the black holes. We also found the entropy expression which does 
not obey the area law for the f(R) black holes. We obtained a Smarr-type formula for the mass, 
and verified that the conserved and thermodynamics quantities satisfy the first law of black hole 
thermodynamics. We also studied the phase behavior of the higher dimensional charged f(R) black 
holes. We found that the system is always thermally stable in four dimensions, while in higher 
dimensions, there is a phase transition in the presence of the conformally invariant Maxwell field 
in R + f(R) theory with constant curvature. 

Finally, we would like to mention that the n-dimensional charged f(R) black hole solutions 
obtained here are static. Thus, it would be interesting if one can construct charged rotating black 
brane/hole solutions from R+f(R) theory in the presence of conformally invariant Maxwell source. 
These issues are now under investigation and will be appeared elsewhere. 
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